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Previous measurements of a quadrupole component of azimuth correlations denoted by symbol V2 
have been interpreted to represent elliptic flow, a hydrodynamic phenomenon conjectured to play 
a major role in noncentral nucleus-nucleus collisions. V2 measurements provide the main support 
for conclusions that a “perfect liquid” is formed in heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC). However, conventional V2 methods based on one-dimensional (ID) azimuth 
correlations give inconsistent results and may include a jet contribution. In some cases the data 
trends appear to be inconsistent with hydrodynamic interpretations. In this study we distinguish 
several components of 2D angular correlations and isolate a nonjet (NJ) azimuth quadrupole denoted 
by V2{2D}. We establish systematic variations of the NJ quadrupole on yt, centrality and collision 
energy. We adopt transverse rapidity yt as both a velocity measure and as a logarithmic alternative 
to transverse momentum pt- Based on NJ quadrupole trends we derive a completely factorized 
universal parameterization of quantity V2{2D}{yt, b, ^/snn) which describes the centrality, yt and 
energy dependence. From j/t-differential V2{yt) data we isolate a quadrupole spectrum and infer a 
quadrupole source boost having unexpected properties. NJ quadrupole V2 trends obtained with 2D 
model fits are remarkably simple. The centrality trend appear to be uncorrelated with a sharp 
transition in jet-related structure that may indicate rapid change of Au-Au medium properties. 
The lack of correspondence suggests that the NJ quadrupole may be insensitive to such a medium. 
Several quadrupole trends have interesting implications for hydro interpretations. 

PACS numbers: 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 


I. INTRODUCTION 


Measurements of the quadrupole component of az¬ 
imuth (/) correlations from RHIC heavy ion collisions in 
the form V 2 = (cos(2 0)) relative to estimates of the 
A-A reaction-plane angle are conventionally interpreted 
to represent elliptic flow, a conjectured hydrodynamic re¬ 
sponse to pressure gradients in the initial collision system 
corresponding to the overlap eccentricity of colliding nu¬ 
clei [T] . In a hydrodynamic (hydro) context [5H1] inferred 
large elliptic flow values combined with other measure¬ 
ments are interpreted to imply rapid thermalization and 
production of a QCD medium with large energy density 
and small viscosity described as a “perfect liquid” m- 

However, questions persist concerning V 2 measure¬ 
ments and interpretations. Conventional ID V 2 meth¬ 
ods [8119] may not distinguish accurately between a non¬ 
jet (NJ) azimuth quadrupole {cylindrical multipole uni¬ 
form on pseudorapidity p over a significant interval near 
midrapidity) and certain jet-related angular correlations 
that vary strongly with p near midrapidity [lOIHSj . The 
terms jet-related and nonjet are discussed in Sec. m 

In previous studies we introduced a physical-model- 
independent method to distinguish geometrically between 
a NJ quadrupole and the quadrupole (m = 2) Fourier 
component of jet-related angular correlations dominated 
by a 2D peak centered at the origin on p and (j) dif¬ 
ference variables [niin]. The notation U2{2D} distin¬ 
guishes the quadrupole component derived from model 
fits to 2D histograms from U 2 {method} data inferred with 
conventional ID methods. We observed that pt-integral 


U 2 { 2 D}(&, ytSNN) data follow simple trends described by 
a few parameters over a broad range of centrality and col¬ 
lision energy above 13 GeV. The trends factorize, 

each factor described by a simple function. That analysis 
was complementary to an analysis reported in Ref. |16j 
focusing on jet-related structure. 

In the present study we extend the NJ quadrupole 
program to measurements of j/t-differential 
U2{2D}(yt, 6, yijvjv) also derived from 2D model 
fits. As an alternative to transverse momentum pt we 
introduce transverse rapidity yt as a logarithmic variable 
compatible with relativistic boost measurements. From 
?/t-differential U 2 { 2 D} data it is possible to infer a 
quadrupole source boost distribution common to hadrons 
of several species m I18j . In this study we determine 
the Au-Au centrality dependence of the quadrupole 
source boost. We also infer a corresponding quadrupole 
spectrum common to several hadron species m and 
substantially different from the spectrum for most 
final-state hadrons. Those results offer new insights into 
possible mechanisms for the NJ azimuth quadrupole. 

This paper is arranged as follows: In Sec. [H] we in¬ 
troduce some general correlation analysis methods. In 
Secs. m and m we describe two alternative methods 
for estimating azimuth quadrupole components of angu¬ 
lar correlations. In Sec. [^ we introduce measured yt- 
differential 2D angular autocorrelations (histograms) de¬ 
rived from particle data. In Sec. |VI| we review system¬ 
atic model-parameter trends from model fits to the 2D 
data histograms. In Sec. VH we define the quadrupole 
source boost and determine its centrality dependence. 
In Sec. |VHI| we extract quadrupole spectra and de- 
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scribe the centrality dependence. In Sec. IX we derive 
a universal factorized parametrization of non-jet quan¬ 
tity z;2{2D}(?/t, b, In S ec, we discuss system¬ 

atic uncertainties, and in Sec. |Xl| we present compar¬ 
isons between quadrupole amplitudes derived from 2D 
fits to angular correlations and from other V 2 methods. 
In Secs.|XlI|andpClII|we present discussion and summary. 


iVA, 4>a) [IS]- The 4>a axis is divided into same-side (SS, 
I'/'aI < 7’‘/2) and away-side (AS, 71/2 < |(?!)a| < tt) regions. 

In the present analysis we impose conditions on the 
space {yti,yt 2 ) to establish the full (yt,b, sys- 

tematics of angular correlations, emphasizing the NJ 
quadrupole obtained from 2D model fits via two model 
parameters: per-particle quadrupole amplitude Aq{ 2D} 
or per-pair amplitude i?Q{2D} (terms defined below). 


II. GENERAL ANALYSIS METHODS 

In this study we report measurements of j/t-differential 
V"2^{2D}(?/t, 5) nonjet azimuth power-spectrum elements 
derived from model fits to 2D angular correlations. 
Transverse rapidity yt (defined below) serves as a log¬ 
arithmic measure of transverse momentum pt . 


A. Kinematic measures and spaces 

A-A collisions with impact parameter b produce 
final-state hadrons in cylindrical 3D momentum space 
{pt,ri,(j)), where pt is transverse momentum, rj is pseu¬ 
dorapidity and (j) is azimuth angle. Transverse mass is 
rrit = -\- m\ with hadron mass mt ■ Pseudorapidity 

is defined by ry = — ln[tan(0/2)] {9 is polar angle relative 
to collision axis z), and ry « cos{9) near ?y = 0. Trans¬ 
verse rapidity is defined by yt = ln[(TOt -\- pt)lrnh]- For 
identified hadrons the proper hadron mass is used and 
yt is then a velocity measure appropriate to test flow 
conjectures. For unidentified hadrons yt with pion mass 
assumed (about 80% of hadrons) serves as a logarithmic 
measure of pt, and the pion mass regularizes the loga¬ 
rithmic trend for small values of pt- The STAR TPC 
acceptance pt > 0.15 GeV/c corresponds to yt > 1. 

Two-particle correlations are structures in the pair 
density on 6D momentum space (j/ti, ?yi, </>i, ?/i 2 ,?y 2 , ^ 2 )- 
Angular correlations can be measured on subspace 
{rii,r] 2 , 4 >i^ 4 > 2 ) given some conditions on transverse mo¬ 
mentum (j>ti,Pti) or transverse rapidity {yti,yt 2 )- We 
can integrate over all yt (?/t-integral analysis) or define 
conditions on two-particle {yti,yt 2 ) (y/t-differential anal¬ 
ysis). Alternatively, we can integrate over some part of 
the angular acceptance (angular acceptance conditions) 
to study conditional correlations on {yti,yt 2 ) [IHHlll- 

An autocorrelation on angular subspace {xi,X 2 ) is de¬ 
rived by averaging pair density p{xi,X 2 ) along diagonals 
on {xi,X 2 ) parallel to the sum axis ccs = xi-\-X 2 [23]. The 
averaged pair density p{xa) on defined difference variable 
XA = Xi — X 2 is then an autocorrelation [53]. For cor¬ 
relation structure approximately independent of xs over 
some limited acceptance Aa: (stationarity, typical over 
27r azimuth and within some limited pseudorapidity in¬ 
terval A7y) angular correlations remain undistorted (no 
information is lost in the projection by averaging) |24j . 
Within the STAR TPC acceptance Ary = 2, Arf = 2 tt [25] 
2D angular autocorrelations are lossless projections of 6D 
two-particle momentum space onto angle difference axes 


B. Correlation structure and interpretations 

The 2D data histograms that form the basis for this 
study (within the STAR TPC acceptance) exhibit the 
same three dominant features from p-p to central Au-Au 
collisions: (a) a SS 2D peak centered at the origin 
on {Va,4>a), (b) an AS ID peak centered at 4>a = 
and approximately uniform on rjA, and (c) an azimuth 
quadrupole component uniform on rjA- Those three ele¬ 
ments are distinguished in all cases by model fits to 2D 
angular correlations on rjA and ^a- Component (a) is 
well-described by a SS 2D Gaussian in most cases. With 
increasing A-A centrality the SS peak is elongated on 
r]A- If high-p( (trigger-associated) cuts are applied the 
SS peak may develop non-Gaussian tails on rjA- Compo¬ 
nent (b) is well-described by a single AS dipole term. 

In p-p and more-peripheral Au-Au collisions compo¬ 
nents (a) and (b) represent intrajet and interjet correla¬ 
tions respectively: Their amplitudes scale with the num¬ 
ber of binary N-N collisions Nmu as expected for dijets, 
and their forms are consistent with pQCD jet structure 
predicted by pythia m and HIJING [28]. They both 
retain the same forms and follow scaling {N-N lin¬ 
ear superposition) in Au-Au collisions up to 50% central¬ 
ity [16j . Throughout that centrality interval it is there¬ 
fore appropriate to refer to (a) and (b) as jet-related 
structures. In more-central Au-Au collisions (above a 
sharp transition near 50% centrality [T3]) the SS 2D peak 
becomes increasingly elongated on tja, and the peak am¬ 
plitude increases faster than scaling. However, other 
features of the SS peak remain consistent with jet pro¬ 
duction [53] . The AS dipole amplitude closely follows the 
SS 2D peak amplitude and also remains consistent with 
jet expectations. Nevertheless, other (nonjet) interpreta¬ 
tions have been proposed for those structures [301 [3I]. 

Azimuth quadrupole component (c) appears to be un¬ 
correlated with jet-related components (a) and (b), for 
example exhibiting a smooth centrality dependence with 
no evidence of the sharp transition [UlIS]. In that con¬ 
text it is appropriate to refer to (c) as the NJ quadrupole. 
The NJ quadrupole inferred from 2D model fits (rep¬ 
resented by symbol n 2 { 2 D}) is an isolated quadrupole 
structure uniform on p near mid-rapidity with maxima 
at 0 and tt on azimuth. Its form is then consistent with 
conventional expectations for “elliptic flow” if that physi¬ 
cal mechanism is relevant. The SS 2D peak (a) projected 
onto ID azimuth can be modeled as a narrow Gaussian 
with its own Fourier series representation |32j . The SS 
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peak Fourier terms then contribute to any ID Fourier 
description of all correlation structure combined. Result¬ 
ing series terms Vm then include admixtures of elements 
(a), (b) and (c). In this study we show that multiple 
(physical) contributions to such a ID Fourier series can 
be distinguished. We therefore refer to NJ quadrupole 
(c) as the object of the present study distinct from a 
jet-related quadrupole derived from (a) which may be a 
source of systematic error for ID V 2 analysis |33j . 


C. Joint and marginal distributions on (yti,yt2) 

This analysis addresses angular correlation system- 
atics corresponding to various cut conditions on space 
{yti,yt2)- Fig. 0 (left panel) shows the relation be¬ 
tween transverse momentum pt and transverse rapidity 
yt- yt ~ 1 (pt ~ 0.15 GeV/c) is the typical lower pt 
bound defining the TPC acceptance for spectra and cor¬ 
relations. The dotted line (m^/2) exp(yt) « pt demon¬ 
strates that yt accurately represents log(p() over the TPC 
Pt acceptance. The close correspondence down to the pt 
acceptance limit arises because the pion mass is « 0.1 
GeV/c. The grid illustrates the cut system for this anal¬ 
ysis, nine bins on yt with fixed width 5yt = 0.4. The pt 
interval covered by the analysis is [0.16,7] GeV/c. 




FIG. 1: Left: The relation between transverse momentum 
Pt and transverse rapidity yt{T^) (assuming a pion mass). The 
grid shows the uniform yt bin system used for the present 
analysis. The dotted line provides a log(pt) reference for 
comparison. Right: The symmetrized two-particle space 
(yt\,yt 2 )- An element of the joint distribution on (yt\,yt 2 ) 
is shown by the bold squares. An element of the marginal 
distribution on yt is shown by the bold rectangles. 


Figure (right panel) shows the binning on {yti,yt 2 ) 
for the yj-differential analysis. The binning system is 
symmetric about the diagonal. The bold squares illus¬ 
trate an element of the joint distribution on (yti,j/t 2 ). 
By integrating over one axis we obtain the ID projec¬ 
tion or marginal distribution on yt represented by the 
bold rectangles. The marginal format is the basis for the 
yt-differential part of this analysis. Since the cut sys¬ 
tem (equal yt bins) and other aspects of this analysis 
are based on yt we prefer that quantity in the text, with 
occasional references to specific pt values. 


D. Single-particle and correlated-pair measures 

The single-charged-particle (SP) yt spectrum is 
Poiyt,b) = cPnch/ytdyt‘^T^drj (azimuth averaged). The 
yt-integral angular density is po{b) = J dytytPo{yt,b) « 
nch/‘2'!^^V averaged over acceptance Ay. The two- 
component particle-yield parametrization [34j as applied 
in this analysis is po{b) = Ppp{Npart/2)[l -b x{i' — 1)] 
where 2'Kppp and x are 2.5 and 0.095 respectively for 
more-central 200 GeV Au-Au collisions. Glauber-model 
centrality parameters Npart/2 and u are defined below. 

p{pi,P 2 ) is the basic pair density on 6D pair momen¬ 
tum space. The event-averaged pair density psib derived 
from sibling pairs (pairs drawn from single events) in¬ 
cludes the correlation structures to be measured, pmix is 
the density of mixed pairs drawn from different but simi¬ 
lar events, pref denotes a minimally-correlated reference- 
pair density derived from (a) a mixed-pair density or (b) 
a product of SP densities via a factorization assumption. 
On pair subspace (yti,yt 2 ) the factorized joint reference 
is Pref{yti,yt 2 ,b) = po{yti,b)po{yt 2 ,b), the marginal ref¬ 
erence is prefivttb) = po{h)po{yt,b) and the yt-integral 
reference is Pref{b) = Po{b)- 

Formation of autocorrelation histograms on difference 
variables (yA,^A) projected from pair angle subspace 
{Vi, 4 ’itV 2 ,(t> 2 ) Fas been described previously [HI 1^ . 
Pair histograms so formed are approximately uniform on 
rjA and (pA- Small deviations from uniformity represent 
correlations of interest. Pair histograms formed by simple 
projection from ( 771 , 772 ) (not autocorrelations) include a 
triangular pair acceptance on difference variable pA- 

Differential correlation structure is determined by com¬ 
paring a sibling-pair density to a reference-pair density 
in the form of difference Ap = psib — Pref representing a 
correlated-pair density or covariance density. There are 
then two choices for a relative correlation measure: 

(a) Per-particle measure Ap/^p^ef has the form of 
Pearson’s normalized covariance wherein the numera¬ 
tor is a covariance and the denominator is the geomet¬ 
ric mean of marginal variances. In the Poisson limit 
a marginal variance is given by = h (x po- Since 
Pref ~ Po ^ Po it follows that the geometric mean of 
variances is given by ytpj^ cx Uch and the normalized co- 
variance is a per-particle correlation measure [TTl [121 US] ■ 

(b) Per-pair measure Ap/pref decreases trivially with 
system size as l/uch- That trend obscures smaller but 
physically-meaningful variations. The per-pair measure 
also tends to increase trivially as a function of yt because 
the pair ratio includes the SP spectrum in its denomi¬ 
nator. The dominant SP spectrum trend also obscures 
physically-meaningful correlation variations. 

In a practical correlation analysis pair ratio Ap/p^e/ —>■ 
Ap/pmix is first calculated directly to cancel particle- 
pair detector inefficiencies. The per-particle measure 
Ap/y/Pref = ^/PrefAp/pmix IS then obtained, where pref 
is constructed from corrected SP spectra po{yt,b) and 
yields po{b) [TB]. 

In this yt-differential analysis we present per-pair 
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quadrupole measurements based on 2D model fits to an¬ 
gular correlations. The basic measures are the Fourier 
components of Ap( 0 a) emphasizing the quadrupole 
term V^- Per-pair ratio Ap/pmix gives BQ{2D}{yt,b) = 
u|{2D}(j/(, &) directly comparable with published V 2 {pt) 
data. Per-particle measure y^PrefAp/p^ix gives 
Aq{2D}(6) = po(^)^^ 2 { 2 D}( 6 ) exhibiting simple sys¬ 
tematic trends on centrality and collision energy m 
[H]. Since Vi{yt,b) = po{b)po{yt,b){b)v^iyt,b) we 
have po{b)AQ{b) = V^{b) and pQ{b)po{yt,b)BQ{yt,b) = 
V^iytjb), defining per-particle Aqib) and per-pair 
BQ(yt,b) quadrupole measures. Per-particle and per- 
pair measures are thus exactly related. For reasons noted 
above V 2 and per-particle measure Aq are the bases for 
physical interpretations. Per-pair measure Bq from this 
2 /t-differential study does not require corrected SP spec¬ 
tra po{yt,b) and provides direct comparison with pub¬ 
lished V 2 (jPt) data. 


a Glauber Monte Carlo. The hatched band, with posi¬ 
tion inferred from angular correlation data, represents a 
sharp transition in jet-related correlation systematics be¬ 
low which jet correlations follow the N-N binary-collision 
scaling expected for linear superposition of N-N colli¬ 
sions (A-A transparency) and no jet modification [16j . 
Fig. i (right panel) shows iz vs fractional impact param¬ 
eter as 1 — b/bo, where bo « 14.7 fm for Au-Au collisions. 

The grids in Fig. [^indicate the centrality bins for this 
analysis defined as follows: Uncorrected minimum-bias 
event samples are divided into 11 nominal centrality bins: 
nine « 10 % bins from 100 % to 10 %, the last 10 % divided 
into two 5% bins. The corrected centrality of each bin 
as modified by tracking and event vertex inefficiencies is 
determined by a running-integral procedure described in 
Ref. [3S]. Centralities from N-N collisions (y « 1.25) to 
central A-A (6 « 0) are thereby determined to 2%. 


E. A-A centrality measures 


A-A centrality is measured by matching the fractional 
cross section cr/tro for some observed Uch to the fractional 
cross section derived from a Glauber Monte Carlo simu¬ 
lation. Glauber parameters Npart/'^ (participant pairs), 
Nbin {N-N binary collisions) and impact parameter b are 
thereby related to Uch integrated within the TPC accep¬ 
tance 1 77 1 < 1. The fractional impact parameter is defined 
by b/bo = a/ct/cto- Centrality measure z/ = 2 Nbin/Npart 
estimates the mean number of N-N encounters per par¬ 
ticipant nucleon (mean projectile-nucleon path length 
across the collision-partner nucleus). We use the same 
Glauber parameters for all energies as purely geometri¬ 
cal measures (the 200 GeV N-N cross section a^N = 42 
mb is assumed for all cases). 



1 - a/a 


0 



l-b/b„ 


F. A-A initial-state geometry measures 


Some features of the initial-state (IS) geometry of A-A 
collisions may influence collision dynamics. IS azimuth 
structure is conventionally modeled by a Glauber Monte 
Carlo. The participant-nucleon azimuth distribution can 
be described by an autocorrelation function on azimuth 
difference c/a [23j. For non-central A-A collisions the 
autocorrelation includes (a) a few even-m sinusoids dom¬ 
inated by TO = 2 (IS quadrupole) phase-correlated with 
vector impact parameter b (the eccentric A-A overlap re¬ 
gion), (b) a uniform background and (c) a delta-function 
term oc Npart (self pairs) uncorrelated with b. 

The Fourier transform of the IS azimuth autocorre¬ 
lation is a power spectrum represented by eccentricity 
elements with per-pair eccentricity mea¬ 

sures [33] 

^m,MC = Cm,opt + fo^ ^ven (1) 

= Se^ for TO odd. 


represents an eccentricity variance due to event-wise 
b fluctuations. Eccentricities em,opt = 2, 4) represent 
the “elliptical” A-A overlap region for fixed b and smooth 
matter distributions. The corresponding to = 2 optical 
eccentricity for 200 GeV Au-Au is parametrized by m 


^2,opt 


1 


log 


10 


^Nbin 


0.96 




0.78 

■ ( 2 ) 


FIG. 2: Left: Mean participant path length v vs fractional 
cross-section measure 1 — (t/cto. The grid shows the eleven 
centrality bins for this analysis. The hatched region indicates 
the “sharp transition” in jet-related correlation properties re¬ 
ported in Ref. m- Right: Path length n vs fractional impact 
parameter b/bg = yja/ao. For Au-Au collisions 6o ~ 14.7 fm. 

Figure[^(left panel) shows participant path length v vs 
fractional cross section in the form l — a/cjo inferred from 


Point-wise Monte Carlo random sampling generates a 
“white-noise” power spectrum i5e^ oc 1/Npart (approxi¬ 
mately uniform on to) corresponding to the self-pair con¬ 
tribution « NpartSi/tA) in the IS azimuth autocorrela¬ 
tion. For a stochastic process there should be no phase 
relation between noise amplitudes and impact pa¬ 
rameter b. All higher to are present in the IS Monte 
Carlo spectrum and might appear in the observed final 
state to some extent i/Monte Carlo sampling at a; « 1/3 
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were a legitimate model of IS geometry relevant to FS 
hadron production for ry « 0 and x » 0 . 01 . 

Figure shows centrality trends for m = 2, 3 IS 
power-spectrum elements on participant-nucleon num¬ 
ber Npart (left panel) and mean participant pathlength v 
(right panel). Plotted are optical eccentricity € 2 ,opt (solid 
curves), Monte-Carlo eccentricity 62 ,mc (dash-dotted 
curves) and so-called “triangularity” Jea (dashed curves). 
From Eq. Q (and ignoring a possible contribution) 
we have e| + Se^ with ~ A/Npart and 

^3,MC “ ^^3 ~ ^/^part- 




FIG. 3: Centrality trends for optical and Monte Carlo 

Glauber initial-state azimuth multipoles with m = 2, 3, on 
participant-projectile-nucleon number (left panel) and binary 
N-N collisions u per participant-nucleon pair (right panel). 


Whether point-like-nucleon sampling at nucleon mo¬ 
mentum fraction a; « 1/3 represents an IS geometry with 
significant manifestations in FS correlation structure for 
X « O.OI is an open question. The present analysis in¬ 
dicates that optical eccentricity £ 2 ,opt is most compatible 
with NJ t; 2 { 2 D} data obtained from 2D model fits. 


III. NONGRAPHICAL NUMERICAL METHODS 


In this section we present a simplified analysis method 
confined to ID azimuth difference /)a corresponding to 
projection p(?7a,/>a) —> p(/’a) of all 2D angular corre¬ 
lation structure within some detector 77 acceptance Ary. 
This description is directly related to conventional ID V 2 
analysis via nongraphical numerical methods (NGNM). 
In Sec. |l^ we describe a more general analysis method 
based on model fits to 2D data histograms on {r]A,4'A), 
the basis for the present NJ-quadrupole study. 


A. Pair densities on azimuth difference <j)A 


can be described exactly by a Fourier series 


I 

PiM = 2 nAn'^ ^ (3) 

00 

= (5((/>A)po/Ary-b -f 2 ^ U^cos(77r(/A), 

m—1 


with = Vqv'^ and Vq = n(n — l)/(27rA?y)^ (over¬ 
line indicates event-ensemble average). The sum upper 
limit is simplified by Uch n. That expression describes 
all pairs within the angular acceptance, including self 
pairs i = j. The Fourier coefficients constitute the 
power spectrum of the azimuth pair density |36j . Any 
and all significant correlation structure projected onto az¬ 
imuth should be accurately represented by the , which 
are additive or extensive measures (whereas the are 
not). The delta function on the RHS represents self pairs 
{i = j) and has its own Fourier series representation (a 
uniform or “white-noise” power spectrum on m) |3fij . 

In this analysis we consider the azimuth quadrupole (a 
cylindrical multipole) term of the power spectrum. The 
m = 2 Fourier coefficient (also a 4D angular density) 
with self-pairs excluded is determined from data by 


I 

yi{b) = d(l)APi4'A,b)cos{2<j)A) 


(4) 


n,n—1 

(27rA7y)2 

= Pl{b)vl{b) ~ pt{h){cos{2(j)A))- 


^ cos[2{(j}, - cjj)]) 

i^i=l,GA?7 

— „2/ 


The approximation in the last line is for n ^ 1. For con¬ 
ventional notation U 2 {method} Eq. ([^ represents V 2 ^{ 2 } 
derived from two-particle ID azimuth correlations. Re¬ 
cent “higher-harmonic” Fourier analysis extends to m > 
2, interpreting any Vm as representing a flow [30]. Equa¬ 
tion 0 is actually equivalent to a model fit to 2D angular 
correlations projected onto ID azimuth. Model-fit com¬ 


parisons are discussed further in Sec. IV 


B. yt-differential Vi{2}(yt,b) measurement 


Figure (right panel) shows joint and marginal con¬ 
ditions on {yti,yt 2 ) (bold lines). The joint distribution 
( 6 D density) on {yti,yt 2 ) is defined by 


p{4>a) is the (4D angular) pair-density projection onto 
azimuth difference i/a obtained by averaging over ryA 
within acceptance Ary. Any distribution on periodic 4>a 


Vi{yti,yt 2 ,b) = 


‘'Vtl 


cos[2((/)i - (j}j)] 


((5?yt27rAry)2 ^ 


ytiyt2 


,(5) 
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where 5yt = 0.4 is the bin width on yt (uniform widths 
in this analysis). The marginal distribution is defined by 


1 

yiiyt,b) = d(j)Ap{yt,(l)A,b)cos{2(j)A) ( 6 ) 


tItj. ,n— 1 , , 

_ 1 cos[ 2 ((^, - (Jjj)] 

Syt(2TTAT])'^ ^ yt 

POO 

= / dy'ty'tVi{yt,y'ub) 

= V2{b)V2{yt,b) 

= Poib)po{yt,b)vl{yt,h), 

where pQ{yt,b) is the SP spectrum on yt- V 2 (yt,b) = 
V 2 {yt,b)/V 2 {b) with V 2 {b) = yJV^ih) defines a self- 
consistent extensive measure system. 

In some conventional NGNM V 2 analyses the pair ratio 
is calculated directly as the ensemble mean of a ratio 


yUyt,b) 


f Ly ^ ^ I L _L 

n„{Ll) S »s|2(*-0,)| (7) 

V2{b)v2{yt,b), 


whereas in other analyses a ratio of mean values is em¬ 
ployed. The pair ratio does cancel detector imperfections 
to some extent. Quadrupole amplitudes V^iyti b) can be 
approximated by using corrected yt spectra poiyt, b) and 
angular densities po{b) averaged over 27r and some rj ac¬ 
ceptance Ay (e.g. STAR TPC). Data from NGNM ap¬ 
plied to two-particle correlations are denoted by U 2 { 2 }. 
The W 2 { 2 D} method used in the present study and based 
on model fits to 2D angular correlations is described next. 


IV. MODEL FITS TO 2D ( 77 , 0) CORRELATIONS 


though a Fourier series can represent accurately any ID 
azimuth distribution, multiple physical mechanisms may 
then contribute to any single Fourier amplitude (azimuth 
multipole) [37]. ID Fourier analysis may be unsuited to 
describe 2D angular correlations from p-p and peripheral 
A-A collisions which include strong variations on pA ac¬ 
curately described by a combination of ID and 2D peaked 
functions mma. In that case a ID Fourier series can¬ 
not describe angular correlations comprehensively over 
all A-A centralities as required to understand the cen¬ 
trality evolution of collision phenomena |32] . 

A 2D data model composed of elementary functions 
can remove ambiguities arising from ID projection onto 
azimuth. The basic premise is as follows: Within some 
limited t] acceptance 2D structure is separated into what 
is approximately uniform on y and what is strongly vary¬ 
ing. The nearly-uniform components are candidates for 
ID Fourier representation. However, alternative repre¬ 
sentations (e.g., ID Gaussian on azimuth) are also con¬ 
sidered. Any components nonuniform on y should be 
modeled by the simplest combination of elementary func¬ 
tions sufficient to describe the 2D data accurately. 


B. Isolating and modeling 2D correlation structure 

If Eq. @ is evaluated for particle pairs drawn from 
the same event the result is the sibling-pair density psib- 
If pairs are drawn from different but similar events the 
resulting mixed-pair density pmix is approximately equal 
to the factorized reference density p^e/ = ^^re/ — Po- 
that context the density of “correlated pairs” (a covari¬ 
ance density with self pairs exeluded) is 

00 

^pif^A) = Psib - Pref = AVq -1- 2 ^ cos{m(j)A)l8) 

m—1 


In this section we extend NGNM analysis of ID az¬ 
imuth projections to graphical analysis of 2D angular 
autocorrelations via fits to data with a 2D model func¬ 
tion including several elements. The y dependence of 
angular correlations is used to distinguish among sev¬ 
eral functional forms that are later interpreted physically 
by comparisons with theory predictions. A per-particle 
model appropriate for studying Glauber linear superposi¬ 
tion in the context of jet production [T3j with amplitudes 
denoted by quantities Ax was applied in a previous yt- 
integral W 2 { 2 D} analysis [T4l|T5]. A per-pair model with 
amplitudes denoted by Bx is applied in the present yt- 
differential analysis to provide direct comparison with 
V 2 {pt) data from conventional ID analysis. The for¬ 
mat does not require corrected SP spectra po{yt,b). 

A. Motivation for 2D model fits 

NGNM data (e.g., V2{2}) derived from ID projections 
discard the y structure of 2D angular correlations. Al¬ 


and the per-pair measure of correlated pairs is pair ratio 
A °° 

— = Avl + 2'^vl, cos(m^A)- (9) 

Note that vf in that series represents a cylindrical mul¬ 
tipole, not the spherical multipole associated with “di¬ 
rected flow” vi [ 3 . The quantity AFq^ = y^ — is 

proportional to fluctuation measure af — n (variance dif¬ 
ference) representing number angular correlations with 
characteristic lengths comparable to or exceeding the ac¬ 
ceptance scale dll US). We wish to extend the mathemat¬ 
ical representation of Eq. ([^ to 2D angular correlations. 

In all 2D data histograms we observe an AS structure 
(AS ridge) that is broad on azimuth and approximately 
uniform on yx within the TPG angular acceptance. The 
latter property implies that the AS structure is a candi¬ 
date for ID Fourier representation, but the former prop¬ 
erty implies that only a few terms in the Fourier series of 
Eq. § are required by the data, and two terms (dipole 
and quadrupole) are sufficient in most cases. An m = I 









7 


AS dipole is generally consistent with pQCD jet structure 
in minimum-bias angular correlations, and the m = 2 NJ 
quadrupole is the object of the present study. 

The remaining 2D structure is strongly varying on 77 A 
and therefore not suitable for ID Fourier series represen¬ 
tation. The 77 -dependent 2D structure is represented by 
a non-Fourier (NF) term. The model function is then 

2 

^piVA,4'A) = Apnf(??a>a)+2 E V^cos{mcj)A)XlO) 

m—1 

where ApNF is a combination of ID (on 77 a) and 2D 
peaked functions plus constant offset. 

In p-p collisions the NF contribution dominates 2D an¬ 
gular correlations and consists of elements predominantly 
associated with either like-sign charge pairs or unlike-sign 
pairs (whereas the NJ quadrupole is observed to include 
both types equally) [^[55]. NF includes a ID peak on 
77 A nearly uniform on (pA and a complex SS 2D peaked 
structure at the {tja, (Pa) origin. Most of the SS 2D peak 
can be modeled by a 2D Gaussian consistent in its prop¬ 
erties with pQCD jet expectations [55]. A smaller 2D 
exponential contribution is consistent with Bose-Einstein 
correlations (BEC) plus conversion-electron pairs [16]. 

In more-central Au-Au collisions the minimum-bias SS 
2D peak becomes elongated on riA and slightly narrower 
on (pA compared to p-p collisions but remains statistically 
consistent with a 2D Gaussian [1^. A NJ quadrupole 
is visually obvious in more-central Au-Au data [TnidS], 
but the quadrupole component remains statistically sig¬ 
nificant for all Au-Au centralities down to N-N colli¬ 
sions [15j . Thus, parametric evolution of a single 2D 
model function with a few simple elements accurately 
represents all collision systems from p-p to central Au-Au. 
This study presents a description of model properties in¬ 
ferred from 2D histograms by T/i-differential model fits. 


C. yt-differential model function 

The 7 /(-integral 7 ) 2 { 2 D}(&) analysis described in 
Refs. |I41115] employed an ll-parameter model function 
as described in Ref. m- For the present j/t-differential 
analysis we introduce three changes to the fit model: 
(a) We switch from per-particle measure Apjy/pref to 
per-pair measure Ap/pref to facilitate comparisons with 
published V 2 {pt) data from ID methods while not requir¬ 
ing corrected SP spectra po{yt,h). (b) We eliminate one 
NF model element not required to describe more-central 
Au-Au collisions, (c) We introduce a Gaussian alterna¬ 
tive model for the AS ID peak. 

With yt cuts imposed jet-related peak structures may 
become narrower, as expected for jet correlations and 
observed in so-called trigger-associated dihadron correla¬ 
tions [38l l39] . If the AS ID peak on azimuth narrows 
its Fourier series representation may require more than 
a single AS dipole term, and an AS ID Gaussian may 
be a more efficient representation. A narrower AS ID 


peak may also introduce a systematic ambiguity between 
the NJ quadrupole and a jet-related quadrupole contri¬ 
bution from the narrower AS peak [JD]- In this analysis 
each of two AS ID peak models is included alternately in 
2D model fits. Any differences in inferred i 3 Q{ 2 D}( 7 /t, b) 
values provide an estimate of systematic uncertainties. 

Given that context we simplify the 2D data model for 
several reasons: (a) The j/t-differential analysis includes 
99 histograms each for 62 and 200 GeV and for two AS 
peak models. Each of nearly 400 conditions then requires 
up to 1000 fits with random starting parameters to insure 
location of global minima. The entire analysis program 
is executed several times to investigate data quality, al¬ 
ternative data models and overall fit stability. Of order 
one million fits are then required, (b) With data subdi¬ 
vided into 11 centrality bins and 9 yt bins the statistical 
power for each fit is reduced, (c) In this analysis we em¬ 
phasize quadrupole systematics for more-central Au-Au 
collisions where some of the NF terms in Eq. (10 1 are not 


required for accurate data description. For those reasons 
a model with fewer parameters is permitted and provides 
improved fit stability and more rapid convergence. 

The per-particle model in Ref. m includes eleven pa¬ 
rameters, of which only five or six parameters represent 
physically relevant model elements [SS 2D peak (3), AS 
ID peak (1 or 2), NJ quadrupole (1)]. The remainder are 
mainly required for structure prominent only in periph¬ 
eral A-A and p-p collisions. The soft-component term 
Ago/i exp{— 77 ^/ 2 CTs} (two parameters), included in the 
present study for more-peripheral collisions, is not re¬ 
quired for more-central Au-Au collisions. In T/j-integral 
analysis the amplitude of that term is observed to fall 
to zero above mid centrality (50% fractional cross sec¬ 
tion) [TB]. We drop the model element representing BEG 
-I- electron pairs (three parameters) because that struc¬ 
ture becomes very narrow in more-central Au-Au col¬ 
lisions. To compensate, three histogram bins near the 
{Fa-:(Pa) origin are removed from the 2D fits to elimi¬ 
nate sensitivity to the narrow BEG -I- electron-pair com¬ 
ponent. 

The 2D model function for j/j-differential analysis ap¬ 
plicable to more-central A-A collisions is then m 


Ap 

Pref 


= Bn 


B2D expi 


(pA 


PA 




-\- Bnll-i- cos(^A - 7r)}/2 -|-Rq 2 cos(2 (/)a),(I1) 

where Bq is a constant offset and Bq —>■ Bq{2D} = 
7 ;|{ 2 D} denotes the NJ quadrupole derived from model 
fits to 2D angular correlations. The soft-component 
term with amplitude Rsoft does not appear explicitly in 
Eq. @ but is included in model fits to more-peripheral 
data histograms. Equation (111 is then referred to as the 
eight-parameter model. 

As noted, the AS ID peak can be modeled by an AS 
ID Gaussian or by its limiting case, an AS dipole [as 
in Eq. If the AS dipole is chosen the in¬ 

ferred quadrupole amplitude determines an upper limit 








on the true nonjet quadrupole. If the AS ID Gaus¬ 
sian is chosen the inferred quadrupole determines a lower 
limit. The difference estimates systematic uncertainties 
for z;|{2D}(j/(,5). The 2D fit model also includes image 
Gaussians at 27r for the SS 2D Gaussian and at —tt for 
the AS ID Gaussian if that peak model is employed [40] , 
Amplitudes are denoted by symbol Bx in the 2D fit 
model of Eq. (11) applied to per-pair data histograms 
Ap/pref as opposed to amplitudes Ax for per-particle 
histograms. From model fits to 2D angular correlations 
we measure pair ratio BQ{2D}{yt,b) = u|{2D}(?/t, 6 ) 
and may convert those data to power-spectrum elements 
via corrected SP yt spectra po{yt,b) (3D densi¬ 
ties) and yields po{b) (2D densities). 


V. 2D ANGULAR AUTOCORRELATIONS 


In this section we introduce data histograms as 2D 
angular autocorrelations. The basic analysis procedure 
is described in Ref. m- Gharged hadrons from Au-Au 
collisions at ^/sxn = 62 and 200 GeV accepted for 
this analysis fell within a detector acceptance defined by 
Pt > 0.15 GeV/c, I 77 I < 1.0 and 27r azimuth. Charge signs 
were determined but particle identification was not oth¬ 
erwise implemented. Further details of track definitions, 
efficiencies and quality cuts are described in Ref. [16]. 
Data for this analysis were selected from earlier RHIC 
running periods where low luminosities insure reduced 
pileup distortions in 2D angular correlations. Reduction 
of pileup effects and other tracking details are described 
in Ref. [TBj. In this presentation we emphasize the 200 
GeV data; the results for 62 GeV are similar modulo the 
energy-dependence factor reported in Refs. [CTIT5] . 



FIG. 4: (Color online) 2D data histograms from 40-50% 
central 200 GeV Au-Au collisions for four yt bins: (a) yt < 1.4, 
(b) 1.8 <yt< 2.2, (c) 2.6 < yt < 3.0, (d) 3.4 < yt < 3.8. The 
narrow BEG -I- electron peak at the origin decreases to zero 
amplitude for larger yt- 


corresponding element. Nonjet quadrupole data are in¬ 
sensitive to that issue and are retained for all yt bins. 


B. Example 2D model fits to yt-difFerential data 


A. yt-difFerential 2D histograms 

Figure [^ shows example yj-differential 2D angular 
autocorrelations for nominal 40-50%-central 200 GeV 
Au-Au collisions. The mid-central case is chosen to il¬ 
lustrate typical model fits over a range of yt bins. The 
?/t-differential analysis is based on 6.7M and 14.5M (both 
year 2004) Au-Au collisions at t/sxN = 62 and 200 
GeV respectively. The panels represent yt intervals (a) 
yt < 1.4, (b) 1.8 < yt < 2.2, (c) 2.6 < yt < 3.0 and (d) 
3.4 <yt < 3.8 {yt bins 1, 3, 5, 7). 

The BEC-electron peak dominates correlations in the 
lowest yt bin as expected. That peak is wide enough in 
the lowest bin that it interferes with the SS 2D jet-like 
peak and degrades the fit quality. Since the information 
obtained on jet structure in the lowest yt bin is mini¬ 
mal, fit values for SS 2D peak parameters from that bin 
are omitted from the rest of the analysis. For larger-?/( 
bins the BEC-electron peak is reduced in amplitude and 
widths. A few histogram bins nearest the origin are re¬ 
moved from the fits and the 2D fit model includes no 


Figure compares fits (left panels) with data his¬ 
tograms (right panels) for two yt bins from 40-50% cen¬ 
tral 200 GeV Au-Au collisions. That figure presents two 
of 99 cases (11 centralities x 9 yt bins). For each case 
approximately 1000 fits starting with randomly-chosen 
initial parameter values are performed. The fit corre¬ 
sponding to the global-minimum value is then cho¬ 
sen. The fits were performed with the eight-parameter 
model function of Eq. 0 including an AS dipole ID 
peak model or with a nine-parameter model including an 
AS Gaussian. The three bins at the origin containing the 
BEG -I- electron peak were excluded from all fits. That 
contribution is mainly confined to the lowest yt bin. 

Small irregularities appearing near the pair-acceptance 
boundary |ryA| = 2 have little effect on 2D fits because 
the statistical uncertainties are largest there. The SS 2D 
peak may deviate from an ideal 2D Gaussian for some 
applied pt cuts, but for this analysis of the marginal dis¬ 
tribution on yt the deviations are not substantial. The 
NJ-quadrupole component is found to be insensitive to 
such deviations (see Sec. for further discussion). 
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FIG. 5: (Color online) Fits to 2D data histograms from 40- 
50% central 200 GeV An-An collisions for two yt bins: (a), 
(b) 1.8 < yt < 2.2 and (c), (d) 3.4 < yt < 3.8. The model 
fits appear on the left in each case. The vertical scale is the 
same for fits and data. The upper limit has been adjusted to 
reveal the BEG -I- electron contribution excluded from the fit 
(narrow peak in three central bins). 



y, y, 


FIG. 6: Per-pair fit parameters derived from 2D fits with 
Eq. 0 to j/t-differential data histograms from 200 GeV 
Au-Au collisions for eight centrality bins (0-70% central), (a) 
SS 2D peak amplitude B 2 D{yt,b), (b, c) SS 2D peak r.m.s. 
widths {yt, b) and {yt, b), (d) Per-pair quadrupole am¬ 
plitude BQ{2T)}{yt,b). Error bars indicate fit uncertainties. 


VI. yt-DIFFERENTIAL Bx{2D} SYSTEMATICS 


Figure shows the yt and centrality evolution of four 
fit parameters from 200 GeV Au-Au collisions describing 
SS 2D peak properties and the NJ quadrupole amplitude. 
The plots show per-pair SS 2D peak amplitude B 2 D, SS 
peak widths ct^/a ^^.d per-pair NJ quadrupole 

amplitude Bq{2T)}. The data for SS 2D peak proper¬ 
ties are not plotted for the first yt bin {pt « 0.2 GeV/c) 
because correlation structure for that centrality is dom¬ 
inated by the narrow peak representing BEG and con¬ 
version electrons [see Fig.[^(a)]. The yt interval for this 
analysis includes almost all jet correlation structure. The 
?/t-differential SS 2D peak parameters are required for the 


V2{2}{yt,b) O V2{2D}{yt,b) comparisons in Sec. XIG 
These results are consistent with Ref. m- 

The B 2 D and Bq data are dominated by trends gener¬ 
ally expected for per-pair amplitudes: (a) decrease with 
increasing centrality and (b) increase with increasing yt, 
both due to the SP po{yt,b) spectrum factor in the de¬ 
nominators of the Bx- Smaller physically-meaningful 
variations are overshadowed by the dominant per-pair 
trends. Just as for Aq{ 2D}(6) the BQ{2D}{yt,b) data 
exhibit remarkable simplicity, but the simplicity is not 
revealed until we present these results in alternative plot¬ 
ting contexts. Note that Bq{2D} values for 0-5% central 
collisions are consistent with zero for almost all yt bins 
(at both collision energies), with small upper limits. 

In Refs, [m [TS] it was reported that j/t-integral per- 


particle NJ quadrupole amplitude Aq{ 2D}(6, y'SATAr) is 
factorizable, the factors having simple functional forms. 
With the 2 /t-differential data from this study we demon¬ 
strate that the yt or pt dependence of the NJ quadrupole 
is also factorizable, leading to a simple quadrupole 
parametrization accurate over a large kinematic space 
as presented in Sec. |IX| 

The sharp transition in jet-related angular structure 
near 50% Au-Au centrality {v « 3) [T^ does not signifi¬ 
cantly alter the pt structure of the SS 2D peak. Similar 
yt trends are observed down to p-p collisions for SS 2D 
peak and quadrupole despite a substantial change in the 
SS 2D peak angular shape. To further explore the NJ 
quadrupole data systematics we introduce the concepts 
of quadrupole source boost and quadrupole spectrum 


VII. QUADRUPOLE SOURCE BOOST 

The broad source-boost distribution reflecting Hubble 
expansion of a bulk medium assumed in conventional hy¬ 
dro descriptions is contrasted with a narrow quadrupole 
source-boost distribution inferred from V 2 {yt) data. 


A. Theoretical context 

According to hydro descriptions elliptic flow is an az¬ 
imuthal modulation of radial flow corresponding to the 
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IS matter eccentricity of non-central A-A collisions, the 
transverse flows arising from large pressure gradients in 
a dense medium uni]. Within the hydro narrative most 
hadrons (especially below 2-3 GeV/c) emerge by “freeze- 
out” from the monolithic flowing medium, and each 
hadron is therefore associated with a particular medium 
speed or relativistic boost corresponding to its freeze- 
out space-time position. Almost all final-state hadrons 
should then reflect the same hadron-source boost distri¬ 
bution. 

As noted in Ref. [3] if all hadrons emerged from a cylin¬ 
drical shell with fixed radial speed (as a limiting case) 
the SP pt spectrum should exhibit a peak at nonzero pt 
and a minimum at zero, reflecting the radial boost of 
the hadron source. On transverse rapidity yt the spec¬ 
trum alteration would be especially simple: the hadron 
spectrum in the stationary lab frame would be the spec¬ 
trum in the moving boost frame shifted to larger yt- For 
V 2 {pt) the consequence would be negative values for pt 
near zero [TS]. But such trends are not expected due 
to “...a more realistic [hadron source] velocity profile, 
[wherein] the peak in transverse-momentum distribution 
disappears.” The “more-realistic” velocity profile is ap¬ 
proximated by that expected for transverse Hubble ex¬ 
pansion of a flowing bulk medium, a broad distribution 
on radial speed Pt extending from zero to some maximum 
value. But so-called “mass ordering” of V 2 {pt) at lower 
Pt should survive as a manifestation of radial flow. 

In Ref. [IH] it was pointed out that the ratio V 2 {pt)/pt 
for several hadron species plotted vs yt with the 
proper mass for each hadron species reveals a com¬ 
mon source boost distribution for identified hadrons 
from a minimum-bias distribution of Au-Au collisions 
(centrality-averaged result). The factor 1/pt emerges 
from a Taylor expansion of the Cooper-Frye expres¬ 
sion [15] for the thermal spectrum from a boosted source. 
Here we consider the centrality dependence of source 
boosts inferred from V 2 { 2 r)}{yt) data for unidentified 
hadrons from 62 and 200 GeV Au-Au collisions. 


B. Quadrupole source-boost centrality dependence 

Figure (left panel) shows the source-boost cen¬ 
trality evolution of yj-differential unidentified-hadron 
Bq{2D} data from Fig. (d) for 200 GeV Au-Au 
collisions. The plotted quantity is unit-normal ratio 
{^/Pt)BQ{2D}{yt,b)/{l/pt)BQ{2D}{b). The ratio for¬ 
mat removes the j/t-integral V 2 {b) centrality dependence 
reported in Refs. [Til US] [see Eq. 0] bringing data 
trends near the spectrum mean pt {yt ~ 1 - 8 ) into align¬ 
ment at unity. The solid curve represents a centrality- 
averaged pion curve from Ref. f or quadrup ole source 
boost Ayto = 0.6 divided by (l/pt)n 2 { 2 D}( 6 ) « 0.1 that 
describes most of the scaled BQ{2T)}{yt, h) data well. 

Figure (right panel) shows the same procedure ap¬ 
plied to data from 62 GeV Au-Au collisions. The results 
are very similar, with some small quantitative differences 



FIG. 7: Left: Per-pair data B2{2T>}{yt,h) from the present 
analysis scaled as described in the text and compared with 
the minimum-bias pion data (solid curve) from Ref. [18] (left 
panel). The vertical dotted lines mark the TPC acceptance 
edge on yt for pions at pt = 0.15 GeV/c. Right: The same 
procedure applied to data from 62 GeV Au-Au collisions. 


noted below. The i/t-integral values U 2 { 2 D}( 6 ) used for 
both plots are derived from Refs. [H [H] for the two 
energies, and see Sec. XHB In both cases (1/pt) = 2 
(GeV/c)“^ and the quadrupole source boost is consis¬ 
tent with Ayto = 0 - 6 . 

The solid and dashed curves in Fig. [7] are defined by 
F{yt) = 9 {b)C{[l - /3t/tanh(?/t)]/(I - Pt)} exp{-pt/P). 
The expression in curly brackets is determined only by 
relativistic kinematics and Ayto US- The other factors 
are derived from data. Factor G is 1.3 and 1.1 respec¬ 
tively for the two energies. Factor g{b) is defined in 
Sec. IX The product g{b)C for three g{b) values (0.55, 


0.75, 1.0) is noted next to the curves (left panel). Ex¬ 
ponential constant P is 4 ± 0.2 GeV/c for 200 GeV and 
5 ± 0.2 GeV/c for 62 GeV. The increase in P for the 
lower energy may result from the softer SP spectrum in 
the V 2 denominator: the spectrum hard component at 62 
GeV is 60% of the hard component at 200 GeV [HI [29] , 
tending to elevate the plotted V 2 ratio at larger yt- 

The plotting format in Fig. [^ includes prefactor 
1 /pt (lab) derived from the measured lab pt ■ Motivation 
for that factor relates to interest in the yt spectrum in the 
boost frame. In the function F{yt) the kinematic factor 
in curly brackets represents the ratio pt(boost)/pt(lab) 
relating pt in the lab frame to pt in the boost frame as 
derived in Ref. m, Eq. (15). If the kinematic factor were 
removed [i.e. if l/p* (boost) were applied as the prefactor] 
most of the data would follow the dash-dotted curves. 

Within the conventional hydro narrative one should 
expect increasing source boosts in more-central A-A col¬ 
lisions as IS particle and energy densities, and therefore 
density gradients, increase. The data in Fig. [^ suggest 
that all scaled BQ{yt,b) data below 20%-central Au-Au 
centrality are statistically identical in shape. The ra¬ 
tio data follow a simple exponential form that facili¬ 


tates the universal parametrization described in Sec. IX 


Uniformity across most centralities suggests that the 
quadrupole source boost is approximately independent 
of Au-Au centrality. Those conclusions are consistent 
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with more-recent Lambda V 2 {pt) data |43j for central 
Au-Au collisions that show the same source boost as the 
centrality-averaged data [HI US] . We pursue that possi¬ 
bility with differential study of quadrupole spectra. 


VIII. QUADRUPOLE SPECTRA 


We next consider the centrality dependence of azimuth 
quadrupole spectrum shapes above pt = 0.5 GeV/c. Ra¬ 
tio measure '(; 2 {method}(j/t, 5) includes the SP yt spec¬ 
trum pQ{yt,b) in its denominator. The SP spectrum has 
a strong jet contribution (spectrum hard component) |46j 
that should not relate to hydro models and is then gener¬ 
ally extraneous to the azimuth quadrupole problem. De¬ 
pending on the V 2 method the numerator of V 2 {yt , b) may 
also include significant contributions from jets in the form 
of a “nonflow” bias. To study the quadrupole spectrum 
in isolation we remove the jet contributions from numera¬ 
tor and denominator of V 2 by focusing on NJ quadrupole 
amplitude Vi{2D}{yt,b) = po{b)po{yt,b)BQ{2J:)}{yt,b). 

Based on BQ{yt, b) data described in the previous sub¬ 
section and Ref. m we define a unit-normal ratio 


Q{yt,b) 


{l/pt)V2H2B}iyt,b) 
{l/pt)Vi{2D}{b) 
{^/Pt)Po{yt,b)v2{2D}{yt, b) 
{^/Pt)po{b)v2{2I)}{b) X g{b)' 


where the data from the present analysis are of the form 
BQ{2D}{yt,b) = V2{2I)}{b)v2{2D}{yt,b). We seek the 
centrality dependence of the quadrupole spectrum shape 
represented by Q{yt,b). The ad hoc 0(1) factor g{b) in 
the second line is defined and discussed in Sec. IIXI 

Figure shows Q{yt,b) data for seven centrality bins 
of 200 GeV Au-Au collisions derived from pair ratios 
BQ{yt,b) obtained in the present analysis (Fig. §. 2D 
model fits with the AS dipole peak model are preferred 
because those fits are more stable. The Bq data are 
combined with SP spectra pQ{yt,b) and yields po{b) from 
Ref. US] to compute Q{yt,b). The ratio is undefined 
for the 0-5% bin since V2{2D}{yt, b) for that centrality is 
consistent with zero for yt > 2. We observe that within 
the data uncertainties j/t-differential quadrupole data fol¬ 
low a universal spectrum shape above yt = 2 {pt = 0.5 
GeV/c) represented by unit-normal quadrupole spectrum 
Qo{yt) (dashed curve) previously derived from MB PID 
V 2 data in the form of a boosted Levy distribution |18j . 
Qo is not a fit to data from the present analysis. We 
conclude that quadrupole source boost Aj/to ~ 0.6 for 
unidentified hadrons (mainly pions) is approximately in¬ 
dependent of Au-Au centrality. That result is consistent 
with 0-10% Lambda V 2 {pt) data from Ref. |43|. 

There are currently no accurate parametrizations avail¬ 
able for 62 GeV SP spectra po{yti b). Thus, complete re¬ 
construction of quadrupole spectra for the lower energy 
is not possible. Based on the results in Fig. the same 
quadrupole spectrum universality may persist there. 



FIG. 8: Unit-normal quadrupole spectra for seven centrali¬ 
ties of 200 GeV Au-Au collisions derived from Eq. ( [T^ . The 
spectrum shapes are independent of Au-Au centrality: all 
Q{yt,b) coincide within systematic uncertainties. The com¬ 
mon shape denoted by Qo{yt) is well described by a boosted 
Levy distribution with parameters T 2 and n 2 . |18j . 


In Fig. [ 7 ] we observed that the plotted data for the 
three most-central bins in the interval yt > 2 {pt > 0.5 
GeV/c) are suppressed relative to the trend for other cen¬ 
tralities. In Fig. |8| th ose data have been rescaled by fac¬ 
tors g{b) in Eq. (12). The resulting vertical shifts bring 
the spectra into coincidence above yt = 2. The agree¬ 
ment of all quadrupole spectrum shapes with common 
form Qoiyt) over that interval is within data uncertain¬ 
ties. The j/j-integral quantity V2{2D}{b) from Ref. [15] 
is consistent with the t> 2 { 2 D}(j/ 4 ,6) values below yt = 2 
for all centralities (because of the SP spectrum shape), 
but a substantially different V 2 {b) describes the V 2 {yt,b) 
centrality trend above yt = 2 for more-central collisions. 

Results from Ref. m and Fig.|^taken together suggest 
that all quadrupole spectra for any A-A centrality and for 
any hadron species follow universal (5o(w)) in the boost 
frame except for most-central Au-Au collisions where an 
additional reduction factor g{b) < 1 is required. That 
conclusion may be contrasted with arguments for quark- 
number scaling of V 2 {pt) data summarized in Ref. [18j . 


IX. V2(yt,b,^sj;^j:7) PARAMETRIZATION 

Model-parameter trends derived from 2D model 
fits to j/t-differential histograms reveal factorization of 
quadrupole systematics on centrality and hadron yt, the 
factors described by simple functions. 
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Combining above results the j/t-differential quadrupole 
amplitude can be described in factorized form as 


Vi{2Q}{yub) « 

That simplicity becomes apparent only in terms of ex¬ 
tensive measure V^iyt^b) obtained by eliminating the 
SP spectrum in the denominator of intensive measure 
with its strong jet contribution. The simplicity 
of Eq. (131 is unique to the quadrupole spectrum. 


The quadrupole amplitude determined in this study is 


BQ{2^}{yt,b) = vl{2-Q}[yt,b) 

= V2{2V)}{b)v2{2T)}{yt,b), 


(14) 


and we use the U 2 { 2 D}( 6 ) parametrization from Refs. m 
ITK] (Sec. |XII^ to infer U2{2D}(yt, b). Based on Sec. VIII 
the V 2 { 2 D\{yt, b) data can be represented accurately 
by a simple parametrization, just as for the j/j-integral 
case. From above we have the relation V^iytib) = 
Po{b)po{yt,b)vl{yt,b). Rearranging Eq. Q accordingly 
we obtain the parametrization 


Po{Pt,b)/po{b) 


QojPt) \ 

Poipt,b)/poib) J 
exp(-pt/P) X f(yt,b), (15) 


V2{2D}(pt,b) = pt (l/pt}v 2 {2D} (b) 
Qo(Pt) 


which can be compared with V 2 {pt, b) data. The ratio 
of unit-normal spectra in curly brackets is approximated 
with reasonable accuracy by exp(—pt/P). An 0(1) em¬ 
pirical factor f{yt,b) representing deviations from that 
function for more-central collisions is 

f{yt, b) = l + C{b) [erf(y* - 1.2) - erf(1.8 - 1.2)](16) 
0(6) = 0.12-(:z-3.4)/5-[(^-3.4)/2]^ 

Function f{yt, b) decreases with yt from values exceeding 
1 for 2 /t < 2 to a constant value g[h) < 1 for j/t > 2 . 
The curve crosses through unity near yt = 1.8 (pt « 0.4 
GeV/c). Values g{b) = 0.9, 0.75, 0.55 are inferred from 
the data for centrality bins 20-30% through 5-10% (the 
0-5% Bq data for yt > 2 are consistent with zero). 

Figure shows comparisons between Eq. (151 and 
200 GeV data from the present j/t-differential analy¬ 
sis (points). The solid dots and open circles represent 
U 2 { 2 D}{pt, 6 } data from 2D fits with AS ID Gaussian and 
AS dipole peak models respectively. The parametriza¬ 
tion of Eq. ( [I5| (solid curves) describes the {2D} data 
accurately over a large kinematic range. Significant dif¬ 
ferences arising from the choice of AS ID peak model 
appear only for large yt values as expected. Thus, the 
algebraic relations in Refs. Oils] and Eq. (f^ are con¬ 
firmed by comparisons with these U 2 { 2 D}(pt 7 ^ data. 

The dotted curves exclude the factor f{yt,b) and 
therefore have the simple form oc ptexp(—p(/P) with 
P = 4 GeV/c for 200 GeV. The more-peripheral data 
(possibly down to N-N collisions) are consistent with 
that parametrization (follow the dotted curves). NJ 



FIG. 9: Pt-differential W2{2D}(pt,6) trends from 200 GeV 
Au-Au collisions (points) compared to a parametrization of 
i!2{2D}(pt, 6) data given by Eq. ( [l^ (solid curves) for (a) 40- 
50%, (b) 30-40%, (c) 20-30%, (d) 10-20%, (e) 5-10%, (f) 0-5% 
centralities. Solid points are for the AS ID Gaussian peak 
mod el a nd open circles are for t he A S dipole peak model in 
Eq. (111. Dotted curves are Eq. (15 I without factor f{yt,b). 


quadrupole data for three more-central bins, and for 
yt > 2 {pt > 0.5 GeV/c), fall increasingly below the 
trend predicted by the parametrization of U 2 { 2 D}( 6 ) from 
Ref. [H] , and quadrupole data for 0-5% central collisions 
are consistent with zero throughout that interval. 


X. SYSTEMATIC UNCERTAINTIES 

Statistical and systematic uncertainties are discussed 
for pt-differential Bq data and for inferred quadrupole 
spectrum trends. Eor this differential study the choice of 
fit model is a compromise between minimizing systematic 
errors and employing the same model to cover large kine¬ 
matic intervals on yt and centrality v. For this discussion 
we refer to the model elements in Eq. 


A. 2D fit-model elements 

At lower yt the 2D exponential (BEG + electrons) peak 
narrows with increasing centrality while the SS 2D peak 
broadens on rj, insuring accurate distinctions. At higher 
yt the exponential peak amplitude drops rapidly to zero. 
Thus, except for peripheral collisions that element can 
be dropped from the 2D fit model. To minimize system¬ 
atic fit errors from that source the central three bins are 
omitted from all 2D fits (bin errors greatly increased). 

At small ly (peripheral collisions, below the sharp tran¬ 
sition) the SS 2D peak is narrow on ry and (j) and accu¬ 
rately described by a 2D Gaussian for all yt bins. In 
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more-central Au-Au collisions the SS 2D peak broadens 
on r] but remains narrow on (p. For most yt bins the 2D 
peak is still accurately described by a 2D Gaussian. At 
larger yt the peak is distorted on ?7 a, developing non- 
Gaussian tails as observed in trigger-associated analy¬ 
sis [38l [39]. The 2D Gaussian model, while no longer 
bin-wise accurate, does estimate the peak amplitude and 
r.m.s. widths satisfactorily for the present study focusing 
on the azimuth quadrupole component. 

In more-peripheral collisions (and in the j/t-integral 
analysis) the AS ID peak is broad enough to be modeled 
accurately by a single dipole term. In more-peripheral 
collisions and at larger yt the AS ID peak appears to 
narrow. If the width of the AS peak becomes substan¬ 
tially smaller than 7r/2 and the AS dipole model is em¬ 
ployed the quadrupole component of the AS ID peak 
might appear as a bias in the inferred quadrupole ampli¬ 
tude. That bias source can be investigated by replacing 
the AS dipole by an AS ID Gaussian model and refitting 
the 2D data. Any differences in inferred i3Q{2D}(?/t, 6) 
establish the systematic uncertainty from that source in 
the inferred NJ azimuth quadrupole amplitude, typically 
at the few-percent level as with other fit uncertainties. 


B. 2D fit-quality systematics 

Figure shows typical fit residuals for lower- and 
higher-yt bins using the eight-parameter fit model. For 
lower-j/t bins (left panel) the residuals do not contain sig¬ 
nificant large-scale structure and are generally consistent 
with statistical uncertainties. The « 0.8 — 2.5. 

The BEG -I- electron peak at the origin appearing in the 
residuals (not described by the fit model) does not con¬ 
tribute to the because those three bins are excluded 
from the fit (assigned artificially large errors). In the 
higher-yt bin the BEG -I- electron peak is negligible, con¬ 
sistent with the pt dependence of both mechanisms. 

To summarize systematic uncertainty trends for the yt- 
differential analysis we identify three zones on kinematic 
space {yt-, v). Zone A is more-peripheral collisions {v < 2, 
1 —(t/(To < 0.3) for all yt, zone B is more-central collisions 
(iz > 2) for lower yt, and zone G is more-central collisions 
for higher yt- The yt boundary between zones B and G 
is approximately yt = 3.8 {pt ~ 3 GeV/c). In zone A the 
eleven-parameter model function for yj-integral analysis 
from Refs. [HI[16] is required for satisfactory fits. In zone 
B either the eleven- or eight-parameter model function 
provides similar fit quality. 

In zone C the eleven-parameter model is excluded be¬ 
cause of fitting ambiguities between part of the SS 2D 
peak structure and the 2D exponential model element, 
but the simpler eight-parameter model function provides 
an adequate description when three bins near the origin 
are excluded from the fit. As noted, the simpler model 
function is therefore utilized to cover zones B and G re¬ 
ported in this -differential analysis, and zone A is not 
reported. Systematic uncertainties in zone B are small 



FIG. 10: (Color online) Residuals from 2D fits to yt- 

differential histograms from 40-50%-central 200 GeV Au-Au 
collisions for Left: yt G [1.8, 2.2] {pt G [0.4,0.65] GeV/c) and 
Right: yt G [3.4, 3.8] {pt G [2,3.3] GeV/c). The vertical scales 
are half the range of the corresponding panels in Fig. [^ 


and consistent with statistical and fit errors. Uncertain¬ 
ties in G may be significant and are explicitly estimated. 

Figure 10 (right panel) shows small but significant 
residuals structure resulting from the non-Gaussian 
shape of the SS 2D peak appearing in laiger-yt bins (zone 
G): excesses at the origin and near the acceptance bound¬ 
aries on 77 a for ())a « 0. A small excess in the inferred 
quadrupole Bq is also observed (e.g., depression near 
4>a = Ti") due to the non-Gaussian SS peak shape. The 
peak-peak residual quadrupole amplitude is about 0.001 
for Bq « 0.014 [Fig. |6 (d)]. Thus, the relative uncer¬ 
tainty is « 0.0005/0.014 < 5%, comparable to the typi¬ 
cal fit uncertainties there. The inferred quadrupole data 
are generally stable against minor changes in jet-related 
fit model elements or the SS 2D peak shape, except for 
the largest yt values where a substantial systematic un¬ 
certainty (20%) must be assigned to the Bq data. What 
matters most for the NJ quadrupole is simply the pres¬ 
ence in the fit model of a SS 2D peak narrow on azimuth. 


C. Quadrupole spectrum uncertainties 

The SP spectrum parametrization from Ref. [46] is not 
constrained by data below yt = 2. Thus, some systematic 
deviations from the Qo{yt) reference may be due to inac¬ 
curacies in the modeled SP spectrum structure. Also, we 
have approximated the hadron spectrum for this study 
by the pion spectrum alone. Protons and kaons do play a 
significant role in the hadron spectrum shape, and those 
spectra are substantially different from the pion spec¬ 
trum, introducing a further source of systematic bias. 

Figure shows relative deviations of unit-normal 
quadrupole spectrum data Q{yt,b) from universal spec¬ 
trum reference Qo{yt)- The symbols are defined as in 
Fig. i The more-central data are scaled by factor 
g{b) < 1 as in Fig. so that spectrum shapes for yt > 2 
{pt > 0.5 GeV/c) can be compared. The more-central 
data then exhibit large deviations below yt = 2 due to the 
rescaling. Such deviations could represent variations in 
SP spectra below yt = 2, but they may also reflect signifi- 
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FIG. 11: Data-model deviations Q{yt, b) — Qoivt) are shown 
relative to fixed spectrum model Qo(yt)- The upper hatched 
band indicates that the r.m.s. data deviation is less than 5% 
except near the acceptance endpoints. There are systematic 
deviations below 0.5 GeV/c and larger fitting uncertainties 
for the largest yt bin. The lower hatched band provides an 
estimate of systematic uncertainties above yt = 1.8 and the 
range of variations corresponding to V 2 {b) vs V 2 {b) below yt = 
1.8. The dashed curves are explained in the text. 


cant changes in the actual quadrupole spectra with phys¬ 
ical implications. The deviations above j/j = 4 are consis¬ 
tent with fit instabilities and sparse statistics. However, 
we also expect an excess at larger yt because those data 
are derived using the AS dipole model. If the AS ID peak 
narrows at larger yt (expected) the quadrupole compo¬ 
nent of the AS ID peak may then contribute a positive 
bias to the inferred W 2 { 2 D} data. 

We can establish an upper limit on possible source 
boost variations with a Taylor expansion of Q{yt,b) 
about Qoiut) relative to variations in boost Ay to- The 
two dashed curves represent ASyto dlog[Qo{yt)]/dyt, with 
dyto = 0 -0‘2Aytg = 0.012 or 2% of the mean source 
boost. The data lie well within those limits for yt > 1.8 
{pt > 0.4 GeV/c). Boost variations with centrality are 
comparable to or smaller than those observed small de¬ 
viations on yt- The comparison suggests that the mean 
quadrupole source boost does not change by more than a 
few percent over a broad centrality interval. In general, 
the Qo{yt) universal quadrupole spectrum represents the 
Q{yt, b) data well over pt in 0.35-4 GeV/c {yt in 1.8-4). 


XI. COMPARISONS AMONG V2 METHODS 


Substantial differences appear between the NJ azimuth 
quadrupole z;2{2D} derived from model fits to 2D angular 
correlations and U 2 { method} data derived from conven¬ 
tional NGNM (equivalent to ID model fits to azimuth 
correlations) [HI [17] . In this section we provide detailed 
comparisons among several methods and consider possi¬ 
ble sources of observed differences. 


A. Algebraic relation between V2{2} and V 2 {EP} 


To establish the relation between conventional event- 
plane {EP} and two-particle correlation {2} methods we 
first note that U 2 {EP} is defined by |5ll^ 


U2{EP} 


^2,obs 

(cOs[ 2(4'2 - ^r)]) ’ 


(17) 


where the numerator is “observed” V 2 ,obs (defined below), 
and the denominator is described as the event-plane res¬ 
olution. To relate V 2 ,obs to U2{2} we define the m = 2 
(azimuth quadrupole) Q vector by [HI SZl EH] 


Qi — 


‘I’kAt] 


Y, u{24>i) = Q2u(2^2), (18) 


2=1,GAt; 


with unit vectors «(</>) and event-plane angle 412 . The 
vector notation is an alternative to that in Ref. m based 
on complex quantities. Q 2 as defined in Eq. (18) is a 2D 
angular density. We then have 


Qi = Q2-Q2 = ^ + V2^{2} 


Po 

27rAr7 


(19) 


as in Eq. @ but with self pairs included in the first term 
on the RHS. It is notable that V'2^{2} and Q 2 differ only 
by the self-pair term. We then have 


V2,obs = (C0S[2((/)- 4'2)]) 


( 20 ) 


\u{ 2 (j)i)-u{ 2 ^ 2 ,i), or 


2=1 


Q2^‘2.obs — 


1 1 


n,n— 1 

— y 

27r Ary n ^ 


u{2(j)i) ■ u{2(j)j) 


= V2{2}V2{2}, 


where Q 2 = Q 2 \/[n — l)/n, and A 2 y indicates that the 
term is excluded from a sum over j. The summation 
condition j ^ i in Eq. (201 (third line) excludes self-pairs 


from that pair sum but not from Q 2 in Eq. (19) (or Q 2 )- 
From Ref. [HI we obtain the “resolution” measure 


cos2[2(4/2-4/,)] = (21) 

(n - 1) 

_ nvl{2} 

1 -I- nu|{2} ’ 

where nw|{2} serves as a statistical figure of merit for 
ratio V 2 analogous to ap^ {a is a nuclear cross section) 
or Lp^ {L is a beam luminosity) for measurements of 
polarization ratio p. The {EPj o {2} relation is then 


V2{2} 


Q 2 In - 1 

V2,obs 

(cOs[ 2(^2 ^r)])r.m.s. 


U2{EP}. 


( 22 ) 
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Small {EP} o {2} differences may arise from covari¬ 
ances corresponding to non-Poisson multiplicity fluctua¬ 
tions. The “event-plane resolution” correction is required 


because invocation of an event-plane estimate via Q 2 im¬ 
plicitly includes a self-pair contribution. Excluding self 
pairs in Eq. (20) does not remove the Q 2 bias. The 
U 2 {EP} estimate does not necessarily relate to an A-A 
reaction plane—it represents all 2D correlation structure 
including jets. We hereafter refer exclusively to u|{2} or 
V' 2 ^{ 2 } except when introducing published U 2 {EP} data. 


B. Algebraic relation between V2{2} and V2{2D} 

The quadrupole power spectrum element V'2^{2} 
(equivalent to pqAq{2}, PqBq{2 }) represents the total 
azimuth quadrupole component for all angular correla¬ 
tions, including both jet-related structures and any non¬ 
jet structure that might be identified with flows. As 
noted in Sec. |IV A] the ?7 a dependence of 2D angular cor¬ 
relations can be employed to separate unique correlation 
components via 2D model fits, as in Refs. [iniiis] and 
the present analysis. Eor almost all collision conditions 
we observe that the AS structure of 2D angular corre¬ 
lations is uniform on 77 A within I 77 I < 1 and completely 
described by a NJ azimuth quadrupole represented by 
Aq{ 2D} and an AS dipole component. 

The only remaining nontrivial structure observed in 
more-central Au-Au collisions is a SS 2D peak (consis¬ 
tent with intrajet correlations). Because the AS dipole 
is orthogonal to all other multipoles the SS 2D peak is 
the only other significant contributor to total quadrupole 
V'2^{2} = w more-central A-A collisions. The SS 

2D peak per-particle quadrupole amplitude (Fourier co¬ 
efficient) is given by [35] 

2Aq{SS}(&) = F2{a^JG{a^^-Ar,)A2D, (23) 

where A 2 D = P 0 B 2 D is the per-particle amplitude of the 
fitted SS 2D peak with r.m.s. widths (tT^^,cr 0 ^), F 2 is 
the m = 2 Fourier component of a unit-amplitude ID 
Gaussian on azimuth with width 

2Fm{(T<j,^) = exp (-mV^^/2) , (24) 

and G < 1 is a calculated 2D ID rj projection factor 
defined in Ref. [35] . We thus obtain the relation 

Aq{2} = Aq{2D} + Aq{SS} (25) 

plus a small contribution from BEG -|- electron pairs 
in more-peripheral collisions. Jet-related quadrupole 
Aq{SS} may be identified with “nonfiow” [3^ 1331133]. 
Nonjet quadrupole Aq{2D} would correspond to ellip¬ 
tic flow if that phenomenon is relevant. We test that 
relation with results from the present and previous 2D 
correlation analysis and published U 2 { method} data in 
the next subsection. 

Strategies have been adopted to reduce nonflow 
(mainly jet contributions) to V 2 by excluding some parts 


of the nominal ( 771 , 772 ) acceptance from NGNM calcu¬ 
lations [35]. For instance, some 77 A interval centered 
at zero may be excluded from projections onto (pA 
by “estimating the reaction plane” with large -?7 detec¬ 
tors [asED]. The motivation is exclusion of jet-related 
structure Aq{SS} from azimuth projections Aq{ 2} based 
on assumptions about the jet fragment distribution on ry. 

Such 77 pair cuts may be less effective at distinguish¬ 
ing jet-related structure from a NJ quadrupole than 2D 
model fits applied within a more-limited 77 acceptance. In 
more-central Au-Au collisions the SS 2D peak is strongly 
elongated and may develop non-Gaussian tails extending 
over a large rjA interval [38]. The effects of 77 -exclusion 
cuts are then quite uncertain and may have little impact 
on jet-related biases in ?; 2 {method} data [35] . 


C. yt-difFerential data comparisons 


Figure 
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shows published u|{EP}(pt,&) data (open 
circles) compared to u|{ 2 D}(pt, b) data from the present 
analysis (solid points or hatched upper limit) and “non¬ 
flow” prediction u|{SS} derived from characteristics 
of the SS 2D peak measured in this analysis. The 
7 ;|{EP}(pt, b) points are obtained by combining V 2 (b) and 
V 2 {pt,b) measurements from Ref. [9]. The hatched re¬ 
gion in the left panel denotes an upper limit on u|{2D} 
[compare with Fig. (f)]. The bold solid curve in 
the right panel is defined (without the factor 100 ) by 
u|{2} = 7;|{2D} -|- 7;|{SS} per Eq. (25 1 . 
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FIG. 12: Left: Comparison of measured 2D angular corre¬ 
lations represented by V 2 {SS}{pt,b) (dash-dotted curve) and 
U 2 { 2 D}(pt, fo) (hatched region, upper limit) with published 
U 2 {EP}(pt, 6 ) data (open circles) for 0-5% central 200 GeV 
Au-Au collisions. Right: Similar comparison for 5-10% cen¬ 
tral Au-Au collisions showing the close correspondence be¬ 
tween the sum 7;|{2D}(pt, 6 )-|-u|{SS}(pt, b) (bold solid curve) 
and published U 2 {EP}(pt, b) data (open circles) as in Eq. ( |25| . 
The light solid curve through the u^lSD} data (solid dots) 
represents the parametrization in Eq. (15l. 


We find that the measured 7 ;|{EP} « ■c|{2} trend is 
predicted by a combination of w|{2D} data and u|{SS} 
representing the m = 2 Fourier component of the SS 2D 
jet peak projected onto ID azimuth. The dash-dotted 
curves u|{SS}(pt, 6 ) derived from SS 2D peak properties 
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inferred from this analysis can be interpreted as the jet 
contribution to u|{2}. We confirm the trend u|{2} = 
u|{SS} + u|{2D} for j/^-differential data based on the 
detailed rj dependence of 2D angular correlations. There 
is no adjustment to accommodate the ^^{EP} data. 

(left panel) shows data for 10-20% cen- 


Figure 13 


tral Au-Au collisions including similar contributions from 
NJ quadrupole u|{2D}(pt,6) and jet-related quadrupole 
u|{SS}(pt, &). The sum (bold solid curve) accurately de¬ 
scribes the published '(; 2 {EP} data. The parametrization 
of Eq. (15) (dotted curve) describes the z;|{2D}(pt,6) 
data over the entire pt acceptance. 
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FIG. 13: Left: Comparison similar to Fig. for 10-20% 
central Au-Au collisions showing near-equal contributions 
from NJ quadrupole V2{2D}{pt,b) and jet-related quadrupole 
U 2 {SS}(pt, 6). The bold solid curve is the sum of {2D} and 
{SS} data. The bold dotted curve is the parametrization 
of V2{2D}{pt,b) from Eq. (|15| combined with a V2{2D}{b) 
value derived from Refs. Ill [l5J (Sec. pOT^ . Right: The 
quadrupole amplitude of the SS 2D peak (dash-dotted curve) 
derived from Eq. ( |23[ ) with data from Fig. [same as Fig. 
(left panel)] compared to a similar calculation using fixed 
SS 2D peak widths (dotted curve and points). The sub¬ 
stantial differences illustrate the importance of accurately- 
measured SS 2D peak properties for understanding jet biases 
in V 2 { 2 } ~ r) 2 {EP} data . 


Figure 13 (right panel) illustrates the importance of 
accurate jet-related 2D correlation m easu rements. The 
dash-dotted curve is u|{SS} from Fig. 12 (left panel) de¬ 


rived from Eqs. ( |23[ ) and ( |24[ ) based on measured SS 2D 
peak characteristics as in Fig. The dotted curve and 
points represent the same computation with the SS 2D 
peak widths held fixed at = 2.5 and = 0.65 
(?/t-integral values for that centrality). This exercise il¬ 
lustrates that accurate description of the u|{EP} data in 
Fig. (left panel) by the dash-dotted curve relies on full 
employment of measured SS 2D peak properties. 

We learn for instance that relative to the correct dash- 
dotted curve the dotted curve assuming fixed SS peak 
widths is too large at lower pt because the SS peak az¬ 
imuth width is substantially larger there, leading to an 
overestimate of Fourier coefficient F 2 {a^^) in Eq. (24) 
by the fixed-width assumption. The dotted curve is too 
large at higher pt because the SS peak ry-width reduc¬ 
tion is not taken into account. The substantial downturn 
in u|{EP}(pt,6) for 0-5% centrality at larger pt is solely 


due to strong narrowing of the SS 2D peak on rjA above 4 
GeV/c toward the p-p value, as in Fig. |^(b). The combi¬ 
nation of measured amplitude and widths of the SS peak 
from Fig. [^accurately describes the nonflow (jet) contri¬ 
bution w|^S} to V 2 {FP} « w|{2} (dash-dotted curve). 


XII. DISCUSSION 

We consider the implications of differential 
V 2 {‘ 2 D}{yt,b) measurements from this analysis for 
conventional V 2 measurements and for interpretations of 
V 2 data in terms of hydrodynamic flows. 


A. Fit models and interpretation of model elements 


We model j/t-differential 2D angular correlations with 
Eq. ( |TT| ) whose elements are motivated only by struc¬ 
tures directly observed in the data, with no a priori 
physical assumptions. We then interpret the elements 
physically by comparison of data systematics with the¬ 
oretical predictions. Based on measured trends for p-p 
and more-peripheral A-A collisions (95-50% fractional 
cross section) we interpret the SS 2D pea k and AS 
ID peak elements as “jet-related” (see Sec. IIB). All 
other elements are then referred to as “nonjet,” includ¬ 
ing the NJ azimuth quadrupole. The same terminology 
is retained in more-central A-A collisions although a jet 
interpretation may be questioned there. The ordered 
sequence—mathematical modeling followed hy physical 
interpretation—is an essential feature. The jet-related 
and non-jet terminology is complementary to flow-related 
and non-flow terms. However, jet phenomenology is well 
established from extensive HEP measurements and QCD 
theory whereas hydrodynamic flows in high-energy nu¬ 
clear collisions (RHIC and LHC as opposed to the Be- 
valac/AGS energy regime) remain a matter of conjecture. 

Jet-related and NJ quadrupole correlation components 
are separately identified. So-called nonflow bias is associ¬ 
ated with the quadrupole (to = 2) Fourier component of 
the jet-related SS 2D peak. Aq{ 2D} or iJQ{2D} results 
are insensitive to the SS peak shape on y as noted in 
the present study. The essential model element is the SS 
peak ID Gaussian factor on azimuth as noted in Ref. EH, 
where ID Fourier fits without an SS peak element are 
strongly rejected by a Bayesian-inference analysis. 


B. yt-integral Aq{2D} systematics 

The Aq{2D}(6, ^/s^n) data from Refs. [21 |2] re¬ 
veal two interesting features: (a) Data from two energies 
are accurately described by the same centrality varia¬ 
tion (defined below), (b) the energy dependence of the 
quadrupole amplitude (in combination with SPS data at 
17 GeV [2]) scales as log(ys]v)v). When plotted on b/bo 
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the Aq{2T)} data reveal a simple variation closely ap¬ 
proximating a Gaussian function centered on b/h^ = 0.5. 

Figure (left panel) summarizes the measured NJ 
quadrupole energy dependence from Bevalac to high¬ 
est RHIC energy. Aq data values at h/ho « 0.5 
minimize the relative effects of jet (nonflow) contri¬ 
butions to Aglmethod}. We observe a major tran¬ 
sition in the energy trend of per-particle measure 
Aq near 13 GeV, suggesting different physical meeh- 
anisms for the measured NJ quadrupole within the 
two energy regimes [53]. Above 13 GeV the function 
R{^JsNN) = log{ys7v]v/13.5 GeV}/log(200/13.5) (solid 
line) describes the energy dependence, with zero inter¬ 
cept at 13.5 ±0.5 GeV. A similar energy dependence was 
observed for {pf) fluctuations/correlations attributed to 
(mini)jets [T5| consistent with jet-related trends observed 
recently at the LHC |52|. The rate of increase with en¬ 
ergy (line slopes) is six times greater at higher energies 
than the Bevalac/AGS trend (dashed curve). 



(GeV) 



R(^SN)NM„(b) 


FIG. 14: Left: Collision energy dependence of V 2 data 

converted to per-particle measure Aq. The solid points are 
from Refs. [H US and Ref. [51]. They follow a logarith¬ 
mic trend proportional to R{^snn) (see text). The Bevalac- 
AGS trend below 13 GeV is proportional to R'{^snn) = 
ln(yGjvjv/3.2 GeV). The open points are taken from Fig. 4 
of Ref. [4^. Right: The azimuth quadrupole amplitude 
Aq{ 2D} divided by optical e^pt plotted vs energy-dependent 
factor R{.^snn) times 200 GeV An-An binary-collision num¬ 
ber Nbin{b). The Aq{ 2D} data are consistent with simple 
proportionality [Eq. ([^] over three decades. 


Figure 14 (right panel) shows (l/e2^opt)AQ{2D}(6) 
vs R{-^SNN)Nbin, where Nbin is the number of N-N 
binary collisions. The vertical-axis variable is moti¬ 
vated to test expectations for the ratio U 2 /e: the trend 
U 2 /e cx S duch/dr], with S the A-A overlap area for more- 
peripheral A-A collisions should transition to a “hydro 
limit” V 2 /e ~ constant in more-central collisions |53j . 
From the present analysis we conclude that for Au-Au 
collisions Aq{2D} data above 13 GeV are described by 

^q{2D}(^, \^snn) = Poib)vl{2I)}{b, ^/sNN) (26) 
= CR{y/SNN)NMnib) ^2,opt (b), 

where coefficient C is defined by lOOGC = 4.5±0.2. Equa¬ 
tion (26) accurately describes measured j/j-integral az¬ 


imuth quadrupole data in heavy ion collisions for all cen¬ 
tralities down to N-N collisions and all energies above 


y'SATjv ~ 13 GeV. It represents factorization of energy 
and centrality dependence for the NJ quadrupole. The 
2D quadrupole data are also consistent with V2^{2D} = 
PoAq{2D} (X NpartNbinel^optib) mi, a trend that, mod¬ 
ulo the IS eccentricity factor, increases much faster than 
the dijet production rate. The non-zero value V 2 ~ 0.02 
from Eq. (26) extrapolated top-p {N-N) collisions agrees 
with a, p-p color-dipole prediction from QCD theory [54j . 


C. yt-differential Aq{2D} systematics 

The BQ{2I)}{yt,b, ^/snn) data from the present study 
reveal two interesting features: (a) the quadrupole source 
boost distribution is independent of Au-Au centrality 
over a large interval (70% to 5%) and (b) the quadrupole 
spectrum shape is independent of centrality over the 
same interval and very different from the SP spectrum 
shape representing most hadrons. The quadrupole spec¬ 
trum is much colder (90 MeV vs 145 MeV) and does not 
change shape above the sharp transition in jet properties 
as does the SP spectrum (e.g. Raa and “jet quenching”). 
Those interesting trends are not apparent from the sys¬ 
tematics of ratio measure V 2 {pt)- Some implications are 
considered in the next subsection. 


D. Physical implications of Aq{ 2D} factorization 


The present study combined with previous yj-integral 
analysis mi US] reveals two factorizations of the NJ 
quadrupole denoted by U 2 { 2 D}: (a) (6, yGjvfv) factor¬ 
ization (j/j-integral case) above 13 GeV and (b) {yt,b) 
factorization at two energies (j/t-differential case). Such 
factorizations become apparent only in terms of extensive 
correlation measure V 2 ^ with accurate distinction be¬ 
tween the NJ quadrupole and other structure, including 
jet-related SS 2D and AS ID peaks. 

In case (a) (j/t-integral case) we can aid interpretation 
by rearranging Eq. (26) to obtain 


2Aq = CR{y^)ic{b)[Npart{b)el,p,ib)]. (27) 

The LHS per-hadron measure of final-state azimuth 
quadrupole Aq{2D} = poV 2 (momentum space) is math¬ 
ematically analogous to the RHS per-participant IS 
quadrupole measure within the square brackets (config¬ 
uration space). The two azimuth correlation measures 
are simply related by the product of participant path 
length v{b) and energy-dependent factor R{^/snn). The 
quadrupole component of the initial A-A overlap source 
depends only on impact parameter b. Thus, the final- 
state NJ quadrupole of produced hadrons (LHS) is simply 
proportional to the IS quadrupole of the collision partic¬ 
ipants (RHS, as determined by y/sjfN,b,A) over a large 
kinematic domain, including N-N (p-p) collisions. 

A plot of Aq vs 1 — b/bo in Ref. [T3| suggests that the 
centrality dependence in Eq. (26) may depend only on 
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the relative impact parameter b/bo independent of colli¬ 
sion energy. The shape of the quadrupole centrality trend 
may not depend on the absolute size of the collision sys¬ 
tem, only on the relative geometry of intersecting spheres 
independent of atomic number A. Further studies with 
lighter nuclei (e.g., Cu-Cu) may test that hypothesis. 

In case (b) (j/^-differential case) we observe that a 
fixed quadrupole spectrum shape is a universal feature 
of Au-Au collisions over most of the centrality range (70- 
5%), and the inferred source boost distribution is narrow 
with fixed mean value. Both results contrast strongly 
with hydro expectations. In the conventional hydro nar¬ 
rative [S] (i) almost all hadrons emerge from a mono¬ 
lithic flowing bulk medium, and (ii) flows are driven by 
pressure gradients corresponding to large IS energy den¬ 
sities in more-central A-A collisions. Item (i) implies 
that quadrupole spectra should be equivalent to SP spec¬ 
tra (and thus cancel in ratio ^ 2 ), that both phenomena 
should reflect a broad source boost distribution corre¬ 
sponding to Hubble expansion of the bulk medium, and 
there should be a close relation with the systematics of 
“jet quenching” in the medium. Item (ii) implies that 
flow magnitudes should increase strongly with A-A cen¬ 
trality, may be negligible in more-peripheral collisions 
and should correspond with jet-quenching systematics. 

In contrast, measured quadrupole spectrum properties 
suggest hadronization from a cold boosted source, pos¬ 
sibly an expanding cylindrical shell. There is no corre¬ 
spondence with the sharp transition in jet properties ob¬ 
served near 50% fractional cross section, suggesting that 
the quadrupole phenomenon is not related to jet forma¬ 
tion through a dense QCD medium. And the narrow 
fixed boost distribution independent of A-A centrality 
appears to be incompatible with a Hubble scenario that 
would describe explosive expansion of a bulk medium, 
the mean boost increasing with A-A centrality |I8) . 

XIII. SUMMARY 

We have obtained azimuth quadrupole component 
V 2 {‘ 2 D} data from transverse rapidity yj-differential cor¬ 
relations for 62 and 200 GeV Au-Au collisions. Ap¬ 
plication of novel analysis methods to 2D angular cor¬ 
relations permits accurate isolation of a nonjet (NJ) 
quadrupole component with simple systematic properties 
on yt, Au-Au centrality and collision energy. 

Conventional f 2 analysis is based on nongraphical nu¬ 
merical methods equivalent to fitting ID azimuth corre¬ 
lations projected from some pseudorapidity 77 acceptance 
with a single cosine function. In the present analysis fits 
with a multi-element fit model are applied to 2D angu¬ 
lar correlations. The fit model is based on identifica¬ 
tion of certain geometric features in the 2D data with¬ 
out assumptions about physical mechanisms. In p-p and 


more-peripheral Au-Au collisions the data features are 
then characterized as jet-related or nonjet by compar¬ 
isons with theory. Those designations are maintained to 
central Au-Au collisions, although some physical inter¬ 
pretations may be questioned in more-central collisions. 

In this analysis we have identified significant “nonflow” 
bias in published U 2 {method}(yt, 6 ) data, the bias de¬ 
rived mainly from a jet-related SS 2D peak. The bias 
is accurately predicted by separately-measured SS peak 
properties. A variety of strategies has been developed 
previously in attempts to reduce the nonflow (jet) bias 
in conventional V 2 data, but the results are inconclusive. 

The systematics of yt-differential V2{2D}{yt,b, y^s^N) 
data from the present study and published yj-integral 
V 2 { 2 D}{b, -^ysNN) data reveal that the quadrupole power- 
spectrum amplitude V 2 {2D}{yt,b, y/sN^) derived from 
those data is fully factorizable. The separate factors on 
rapidity, centrality and energy are represented by sim¬ 
ple functional forms. In terms of per-particle quadrupole 
measure Aq = (Po is the single-particle density) 
the energy dependence is observed to be proportional 
to log(s/so) (v^ ~ 13 GeV) as expected for a QCD 
process. The centrality dependence is essentially Gaus¬ 
sian on relative impact parameter b/bg. The quadrupole 
power-spectrum centrality trend is l 2 ^{ 2 D}( 6 , y^sjvfv) oc 
NpartNbiniopti^)- The Same trends accurately describe 
data from p-p to mid-central Au-Au collisions. A nonzero 
V 2 value for p-p collisions derived by extrapolation is con¬ 
sistent with a theory prediction based on an alternative 
(nonflow) QCD mechanism for the NJ quadrupole. 

From the yt-dependence factor quadrupole spectra can 
be reconstructed, and a quadrupole source boost inferred 
for each collision system. The quadrupole spectrum 
shape is the same for three hadron species and for all col¬ 
lision systems, and the quadrupole source boost (a single 
value) is approximately independent of Au-Au centrality. 

Our results have implications for hydrodynamic inter¬ 
pretations of A-A collisions. The universal quadrupole 
centrality trend can be contrasted with the trends for jet- 
related correlations which exhibit a common sharp tran¬ 
sition within a small centrality interval, from N-N linear 
superposition in more-peripheral Au-Au collisions to a 
substantially different dependence in more-central colli¬ 
sions. In contrast, the trend for V2{2D}/CopUcai remains 
smooth and slowly varying from p-p to more-central 
Au-Au collisions. If jet production responds to forma¬ 
tion of or changes in a dense bulk medium the azimuth 
quadrupole appears unresponsive to such a medium. 

Further implications for hydro models arise from 
quadrupole spectrum results and quadrupole source- 
boost trends. The mean source boost does not vary sig¬ 
nificantly with Au-Au centrality, and the narrow boost 
distribution is inconsistent with Hubble flow of an ex¬ 
panding bulk medium. The NJ quadrupole and ratio 
7 ; 2 { 2 D}/e fall to zero for most-central Au-Au collisions. 
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